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Abstract 

In this work, we are interested in the small time global null controllability for the viscous 
Burgers' equation yt — yxx + yVx = u(t) on the line segment [0, 1]. The second-hand side is a 
scalar control playing a role similar to that of a pressure. We set y(t, 1) = and restrict ourselves 
to using only two controls (namely the interior one u(t) and the boundary one y(t,Q)). In this 
setting, we show that small time global null controllability still holds by taking advantage of both 
hyperbolic and parabolic behaviors of our system. We use the Cole-Hopf transform and Fourier 
series to derive precise estimates for the creation and the dissipation of a boundary layer. 

1 Introduction 

1.1 Description of the system and our main result 

Let T > be a positive time, possibly small. We consider the line segment x G [0, 1] and the following 
one-dimensional viscous Burgers' controlled system: 



(1) 



The scalar controls are u £ -L^(0, T) and v € i/^/''(0, T). The second-hand side control term u(-) plays 
a role somewhat similar to that of a pressure for multi-dimensional fluid systems. Unlike some other 
studies, our control term u(-) depends only on time and not on the space variable. 

For any initial data ?/o S ^^(0, 1) and any fixed controls in the appropriate spaces, it can be shown 
that system ^ has a unique solution in the space X = L'^{{0,T); H^{0, 1)) nC°([0, T]; L2(0, 1)). This 
type of existence result relies on standard a priori estimates and the use of a fixed point theorem. 
Such techniques are described in [33] • One can also use a semi-group method as in |[35J. Our main 
result is the following small time global null controllability theorem for system ([1]): 

Theorem 1. Let T > be any positive time and yo by any initial data in L^(0, 1). Then there exists 
a control pair u € L°°{0,T) and v G _ff^/'*(0,T) such that the solution y £ X to system {ip is null at 
time T. That is to say, y is such that y(T, •) = 0. 



yt + yvx - Vxx 


^u[t) 


in 


(0,T)x (0,1), 


2/(t,0) 


^v{t) 


in 


(0,T), 


2/(t,l) 


= 


in 


(0,T), 


2/(0, x) 




in 


(0,1). 



1.2 An open-problem for Navier-Stokes as a motivation 

As a motivation for our study, let us introduce the following challenging open problem. Take some 
smooth connected bounded domain SI in R2 or ]r3^ 

Consider some open part F of its boundary d^. 
This is the part of the boundary on which our control will act. We consider the following Navier-Stokes 
system: 

Vt- ^y+{y ■ V)?/ = -Vj3 in 
div y = in 
y = on 
y(0, •) =yo(-) m 
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We consider this system as an underdetermined system. Our control will be some appropriate trace 
of a solution on the controlled boundary F. 



Open problem 1. Is system small time globally null controllable? That is to say, for any T > 
and yo in some appropriate space, does there exist a trajectory of system (0) such that y(T, •) = 0? 

Many works have be done in this direction. Generally speaking, one can distinguish two approaches. 
First, one can think of the nonlinear term as a perturbation term and obtain the controllability by 
means of the Laplacian term. For instance, Fabre uses in (TH] a truncation method for the Navier- 
Stokes equation. In [34 1, Lions and Zuazua use Galerkin approximations for various fluid systems. 
Of course, this approach is very efficient for local results. The most recent result concerning local 
controllability for system ^ is the one contained in [22 by Fernandez-Cara, Guerrero, Imanuvilov 
and Fuel. Their proof uses Carleman estimates. 

The other approach goes the other way around. Indeed, in finite dimension, it is known that if 
y — F{y) + Bu where F is quadratic is controllable, then ij — F{y) + Ay + Bu is controllable too 
(see |151 Theorem 3.8]). Likewise, for fluid systems, trying to get a fast controllability result implies 
to work at high Reynolds number (ie. with big fluid velocities, or low viscosity) inside the domain. 
Therefore, inertial forces prevail and the fluid system behaves like its null viscosity hyperbolic limit 
system. In our case, we expect to deduce results for Navier-Stokes from the Euler sytem. For Euler, 
global controllability has been shown in [10 1 by Coron for the 2D case (see also [T^]) and by Glass 
for the 3D case in [5^. Their proofs rely on the return method introduced by Coron in [5] (see 
also [13, Chapter 6]). For Navier-Stokes. things get harder. In [16J, Coron and Fursikov show a global 
controllability result in the case of a 2D manifold without boundary. In [23J, Fursikov and Imanuvilov 
show a global exact controllability result for 3D Navier-Stokes with a control acting on the whole 
boundary (ie. F = d^). 

Other approaches exist. Let us mention for instance the work |2j, where Agrachev and Sarychev 
control Navier-Stokes equations by means of low modes. They use methods of differential geomet- 
ric / Lie algebraic control theory for finite dimensional control systems. 

The main difficulty of Open problem 1 is the behavior of the system near dQ,\V. Indeed, although 
inertial forces prevail inside the domain, viscous forces play a crucial role near the uncontrolled bound- 
ary, and give rise to a boundary layer. An example of such a phenomenon can be found in [llj where 
Coron derives an approximate controllability result and highlights the creation of a boundary residue. 
Hence, the key question is whether one can handle such a boundary layer by means of the control. 

Some authors have tried to study simplified geometries for Open problem 1. In [7J, Chapouly stud- 
ies a Navier-Stokes equation on a rectangle with Navier-slip boundary conditions on the uncontrolled 
part of the boundary. She obtains small time global null controllability. In [S^ and [55], Guerrero, 
Imanuvilov and Fuel prove approximate controllability for a Navier-Stokes system in a square (resp. 
in a cube) where one side (resp. one face) is not controlled and has zero Dirichlet boundary condition. 

Burgers' equation has been extensively used as a toy model to investigate properties of more 
complex systems in a rather simple setting. This equation was introduced in the seminal paper [S] 
by Burgers. Both from a theoretical and a numerical point of view, it already exhibits some key 
behaviors (such as interaction between the non-linearity and the smoothing effect). Therefore, our 
Theorem [1] can be seen as an example for fast global null controllability despite the presence of a 



dQ.\V 




Figure 1: Setting of the Navier-Stokes control problem 
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Dirichlct boundary layer. Moreover, despite the simplicity of Burgers' equation, the analogy between 
systems ([T]) and ^ is quite striking. We can interpret our scalar control u{t) as some one-dimensional 
counterpart of a pressure gradient for 2D or 3D. 

1.3 Previous works concerning Burgers' controllability 

Concerning the controllability of the inviscid Burgers' equation, some works have be carried out. In [3], 
Ancona and Marson describe the set of attainable states in a pointwise way for the Burgers' equation 
on the half- line x > with only one boundary control at x = 0. In ^31j . Horsin describes the set of 
attainable states for a Burgers' equation on a line segment with two boundary controls. Thorough 
studies are also carried out in [T] by Adimurthi et al. In [33], PerroUaz studies the controllability of 
the inviscid Burgers' equation in the context of entropy solutions with the additional control u{-). 

Let us recall known results concerning the controllability of the viscous Burgers' equation. We 
start with some positive results. 

First, Fursikov and Imanuvilov have shown in [25j a small time local controllability result. It 
concerns local controllability in the vicinity of trajectories of system ([T]) and it only requires one 
boundary control (either y(t, 0) or y{t, 1)). Their proof relies on Carleman estimates for the parabolic 
problem obtained by seeing the non-linear term yy^ as a small forcing term. 

Global controllability towards steady states of system ([T|) is possible in large time both with one or 
two boundary controls. Such studies have be carried out by Fursikov and Imanuvilov in [24] for large 
time global controllability towards all steady states, and by Coron in [14J for global null-controllability 
in bounded time (ie. bounded with respect to the initial data). 

When three scalar controls (namely u{t), y{t,0) and y{t, 1)) are used, Chapouly has shown in [B] 
that the system is small time exactly controllable to the trajectories. Her proof relies on the re- 
turn method and on the fact that the corresponding inviscid Burgers' system is small time exactly 
controllable (see |13l Chapter 6] for other examples of this method applied to Euler or Navier-Stokes) . 

Some negative results have also been obtained. 

In the context of only one boundary control y(t, 1), first obstructions where obtained by Diaz 
in [TH|. He gives a restriction for the set of attainable states. Indeed, they must lie under some limit 
state corresponding to an infinite boundary control y(t, 1) = -t-oo. 

Still with only one boundary control, Fernandez-Cara and Guerrero derived an asymptotic of the 
minimal null-controllability time T{r) for initial states of norm lower than r (see [2T]). This shows 
that the system is not small-time controllable. 

Guerrero and Imanuvilov have shown negative results in [57| when two boundary controls y{t, 0) 
and y{t, 1) are used. They prove that neither small time null controllability nor bounded time global 
controllability hold. Hence, controlling the whole boundary does not provide better controllability. 

1.4 Strategy for steering the system towards the null state 

In view of these results, it seems that the pressure-like control u(t) introduced by Chapouly is the key 
to obtaining small time global controllability results. In order to take advantage of both hyperbolic 
and parabolic behaviors of system ([T]), our strategy consists in splitting the motion in three stages: 

Hyperbolic stage: Fast and approximate control towards the null state. During this very short 
stage t e [0,er] where < e ^ 1, the systems behaves like the corresponding hyperbolic one, as the 
viscous term does not have enough time to act. This hyperbolic system is small time null controllable. 
During this first stage, we will use both u{-) and v{-) to try to get close to the null state, except for a 
boundary layer at a:; = 1. 

Passive stage: Waiting. At the end of the first stage, we reach a state whose size is hard to 
estimate due to the presence of a boundary layer. During this stage, we use null controls v(t) — 
u{t) — 0. Regularization properties of the viscous Burgers equation dissipate the boundary layer and 
the size of y{t, •) decreases. We show that it tends to zero in L^(0, 1) when e 0. This is a crucial 
stage as is enables us to get rid of the boundary residue. It seems to be a new idea and could also be 
applied for other boundary layers created when trying to get fast global controllability results. 

Parabolic stage: Local exact controllability in the vicinity of zero. After the two first stages, we 
succeed in getting very close to the null state. The non-linear term becomes very small compared to 
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the viscous one, and the system now behaves Uke a paraboUc one. We use a smaU time local exact 
controllability result to steer the system exactly to zero. Dm'ing this last stage, we only need the 
control v{-). 



Most of the work to be done consists in deriving precise estimates for the creation and the dissi- 
pation of the boundary layer. We will use the Cole-Hopf transform (introduced in [H] and [30]) and 
Fourier series to overcome this difficulty. First, we will investigate the hyperbolic limit system (see 
Section [5]). Then we will derive estimates for the creation of the boundary layer during our hyperbolic 
stage (see Section [3]) and estimates for its dissipation during the passive stage (see Section |4]) . This 
will achieve the proof of a small time global approximate null controllability result for our system ([1]). 
In Section [5l we will explain the parabolic stage and the local exact controllability. 



1.5 A comparison lemma for controlled Burgers' systems 

Throughout our work, we will make an extensive use of the following comparison lemma for our 
Burgers' system, in order to derive precise estimates. When the viscosity is null, this comparison 
principle still holds for entropy solutions (as they are obtained as a limit of low viscosity solutions) . 

Lemma 1. LetT^v > and consider yo^ij^ G i^(0, 1), u,u £ L^{0,T), wo, wo, "^ii '^i: G H^^'^{0,T). 
Assume these data satisfy the following conditions: 

yo ^ yo o,nd u < u and vq < Vq and vi < vi. 

Consider the following system (which is a generalized version of system 

yt + yyx 

(3) 



- i^yxx 


^u{t) 


in 


(0,T)x (0,1), 


y{t,o) 


= vo{t) 


in 


(0,T), 


yit,i) 


= Vi{t) 


in 


(0,T), 


2/(0, x) 


= yoix) 


in 


(0,1). 



Then the associated solutions y,y G X to system are such that: 

y<y on (0,T)x(0,l). 

One can find many comparison results in the litterature (see for instance the book [37] and the 
references therein). However we give the proof of Lemma [T] both for the sake of completeness and 
because with have not found this precise version anywhere. 

Proof. We introduce w — y — y. Thus, w G X is a solution to the system: 

Wt-Wxx ^ [u - u) ~ \{wy + wy)x in (0,r)x(0, 1), 

w{t,Q) ^vo{t)-vo{t) in (0,r), 

w{t,l) ^vi{t)-vi{t) in (0,r), 

w{0,x) ^yo{x)-yo{x) in (0,1). 

We want to study the negative part oiw: 6 = mm{w, 0). Hence, S{t, 0) ~ 6{t, 1) ~ 0. Now we multiply 
the evolution equation by (5 < and integrate by parts for x E [0, 1] to get a i^-energy estimate for S: 

Thus, we can incorporate the first term of the right-hand side in the left-hand side: 

Since y,y e L'^ {{0,T); H^{0, 1)), we have that: 

t^\\y{t,-)+y{t,-)\\l, belongs to L\0,T). 
Hence we can use Gronwall's lemma. Since 5(0, •) = 0, we deduce that 6 = and y < y. □ 
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2 Analysis of the hyperbolic limit system 



2.1 Small time versus small viscosity scaling 

Let us choose some e > 0. We want to study what happens during the time interval [0,eT]. To study 
this very short first stage, we perform the foUowing change of scale. For t G [0, T] and x G [0, 1], let: 



y{t,x) = ey{et,x). 
Hence, y G X is now the solution to the small viscosity system: 

yt + yyx 



(4) 



(5) 



where we performed the following scalings: u{t) — e^u(et), v{t) — ev{et) and ya{x) — eyo{x). This 
scaling is fruitful because it highlights the fact that, when small time scales are considered, the non- 
linear term is the key term. We want to understand the behavior of the limit system when e = 0. 
Therefore, let us consider that -u(-), v{-) and j/o(-) are fixed data, and let e go to zero. 



^Vxx 




in 


(0,T)x (0,1), 




= m 


in 


(0,T), 


y(i,i) 


= 


in 


(0,T), 


y(0,x) 




in 


(0,1), 



2.2 Obtaining the entropy limit 

When one considers the entropy limit e — )■ for system ([5]), it is not possible to keep on enforcing strong 
Dirichlet boundary conditions. A boundary layer appears and it is necessary to weaken the boundary 
conditions. Otherwise, the system would become over-constrained. The pioneer work concerning the 
derivation of such weak boundary conditions is the one by Bardos, Le Roux and Nedelec in [4J. In 
our particular setting, one gets the following system: 

yt + ^iy% =u(t) in (0,T)x(0,l), 

y{t,0) £E{v{t)) in (0,T), 

y{t,l) >0 in (0,T), 

y{0,x) ^yoix) in (0,1), 

where 

E(a) ^ / if a<0, 

^ ' \ ]-oo;-a]U{a} if a > 0. 

Let us explain the physical meaning of the set E{-). On the one hand, when one tries to enforce 
a negative boundary data on the left side, characteristics instantly flow out of the domain, and our 
actions are useless. On the other hand, if we set a positive boundary data, then: either it is satisfied, 
or a greater negative wave overwhelms it. 

Without getting into the details of entropy solutions (for that subject, refer to the definition given 
in [3] or to the book [38j), we will use the following theorem that guarantees that system ^ is 
well-posed. 

Theorem 2 (Bardos, Le Roux and Nedelec in U]). For any initial data yo e BV{0, 1) and any pair 
of controls u G L^{Q,T), v G BV{0,T), system ^ has a unique entropy solution y in the space 

w((o,i)x (o,r)). 



2.3 Small time null controllability 

We are going to show a small time null controllability result for the hyperbolic limit system. However, 
this will not imply small time global controllability since the system is not time reversible. Indeed, 
even though the PDE seems time-reversible, the definition of an entropy solution is not. 

Theorem 3. System (0) is small time globally exactly null controllable. 

Let us start by giving the intuition of the proof. In a first step, we enforce a constant left boundary 
data H > 0. It moves towards the right and overrides the initial data yo{-) provided that the shocks' 
propagation speed is sufficient. Therefore, H is chosen by using the Rankinc-Hugoniot formula. 
Figure [2] shows a simulation of this first step for some smooth initial data yo. At the end of this step, 
we have y{-) = H. During the second step, we use some constant negative u to get back down to the 
null state. 

Now let us give a rigorous proof using the comparison principle. 
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Figure 2: Overriding of an initial data yo{x) by some constant state y{x) = H for system 



Proof. Let yo{x) £ BV{Q, 1) and T > 0. Let us choose H such that: 

i(i^-||yolk~)>|. (7) 

We enforce the foUowing controls: 

_ i H for t e [0,T/2], 

^ " \ 2F(l-^) for te[T/2,T], 




for te[o,r/2], 

for te[T/2,T]. 



From Theorem [21 we know that there exists a unique entropy solution y £ BV{{0,1) x (0,T)) for 
these data. Let us show that y{T/2, •) = H. Therefore, we will easily deduce y{T, •) = 0. 



Let us extend our initial data from [0, 1] to M. Since Theorem [5] guarantees the uniqueness of the 
solution, the restriction to a; G [0, 1] of our global solution will be the unique solution to ©• Therefore 
we consider yo G BV{M.): 

( H for x < 0, 
^0(2:) = < yai^) for < x < 1, (8) 
[0 for 1 < x. 

Let us introduce y the weak entropy solution defined on M x [0, T] associated to this initial data. 
Thanks to Rankine-Hugoniot formula and (O, we know that: 

y{t,x)=H for x < ^ ■ ~ 

Hence, y(T/2, x) — H ioi x E [0, 1], and y{t, 0+) = H. If we want the restriction of ?; to be a solution 
to ([S]), we need to check that y{t, 1") > 0. Let us use the comparison principle for solutions to inviscid 
Burgers' equation. It can be obtained by taking the null viscosity limit in our Lemma [T] Hence 
y{t,x) > w{t,x) where w is the solution associated to the initial data: 

( H for X < 0, 

wo{x)^< -||j/o|ioo for 0<a;<l, (9) 
[0 for 1 < a;. 
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We have two Riemann problems. Near a; = 1, we have a rarefaction wave. Hence x i— ^ w{t,x) is 
continuous near a; = 1 as long as the H shock wave has not reached x = 1. Hence w{t, 1^) = before 
T* = 1/{2H - 2||j/o||cx)), then w{t, 1") = H. This is why w{t, 1") > 0. Thus 1") > w{t, l") > 0. 
The restriction j/[o,i] is the unique solution to ^ and it is equal to H at time t — T/2. □ 

This proof uses the comparison principle for Burgers' equation. Since we consider a 1-D system, 
this is not a problem. However, if we wanted to be able to handle multi-dimensional systems, we 
could use the generalized characteristics method from Dafcrmos (see fl71). This technique has been 
successfully used by PerroUaz in [55] . 



3 Hyperbolic stage and settling of the boundary layer 

Thanks to the analysis of the hyperbolic limit system, we were able to exhibit controls steering the 
system towards the null state from any initial data. Now we want to apply the same strategy to 
the slightly viscous system ([SJ by using very similar controls. However, a boundary layer is going to 
appear. Our goal in this section is to derive bounds for the boundary layer at the end of this stage. 



3.1 Steady states of system ([5]) 



From now on. the viscosity is positive. Hence, since we have a zero Dirichlet boundary condition 
y(l) = 0, we cannot hope to reach a constant state y{x) = H > . However, we expect that we 
can get very close to the corresponding steady state. Let us introduce the following steady state of 
system ([S]): 

h%x) = iJtanh ( ^(1 



2e 



(10) 



Lemma 2. For any H > and any e > 0, h"^ defined by ilU\) is a stationary solution to system 
with controls: u{t) — and v{t) — _ff tanh (||) • 

Proof. The proof is an easy computation. In fact, it is possible to compute explicitly all the steady 
states for system ([5]), at least when u — 0. This is done in [24J with viscosity e — 1. □ 

We have chosen a boundary data v{t) = 7J tanh (||-) for the definition of our steady state /i^, but 
we will use a control v{t) = H for the motion. This technical trick will lighten some computations and 
is relevant since both terms are exponentially close as e goes to zero. However, some proofs require the 
use of the exact steady state corresponding to a boundary data v{t) = H . Therefore, we introduce: 



k^x) = is: tanh ( ^(1 



where K > is given by the implicit relation K ta.uh {K / (2e)) — H. 



(11) 



Lemma 3. For any H > and any £ > 0, fc^ defined by Ul]) is a stationary solution to system (0j 
with controls: u{t) — and v{t) — H . Moreover, we have the estimate: 



IF 



(12) 



Proof. Lemma m gives us that fc^ is a steady state. For the estimate, we write 









< 


Kt&nh 1 — I 


- iJtanh 1 — I 




\2eJ 


\2eJ 



< H 

< 2He-"'^ 



1 — tanh — 



2e 



□ 
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3.2 First step: overriding the initial data 



In order to get close to the steady state h'^ , it is necessary to choose H in such a way that a Rankine- 
Hugoniot type condition is satisfied. Once we get close enough to the steady state, the solution will 
very quickly converge to the steady state. Indeed, the eigenvalues of the linearized system around this 
steady state are real, negative, and of size at least 1/e. This guarantees very quick convergence to the 
steady state. Such a study of the linearized problem around a steady state for the Burgers' equation 
can be found in [32]. We give the following lemma describing the settling of the limit layer. 

Lemma 4. Let T > 0, H > and j/o € Hq{0, 1) be given data. Then for e > small enough, there 
exists a boundary control v <E H^^'^{0,T) such that v(-) < H and such that the solution y £ X to 
system with initial data yo — syo and controls u — and v satisfies: 



ly(r,.)-/^e(-)llL^(0,l) 



-l/2g-ft(HT-2) 



(13) 



h'^ix) 




Figure 3: Example of evolution from an initial data towards a steady state. 

Let us postpone the proof of Lemma|3]for the moment. We start by giving a few remarks concerning 
this statement and its proof. The intuition is to choose a boundary control v(t) = H , just like we have 
done for the hyperbolic case. Moreover, we want to use the Cole-Hopf transform and Fourier series to 
compute explicitly y{T, •). Let us introduce the Cole-Hopf transform: 



Z(t, x) = exp 



1 

'2e 



yit,s)ds 



This leads to the following heat system for the new unknown Z: 



Z(i,0) 
Z{0,x) 



= 1 

= 

= Z^x) 



^yAt,o))z 



on 
on 
on 
on 



(0,T)x (0,1), 
(0,T), 

(o,r), 

(0,1), 



(14) 



where the initial data Z*^ is computed from the initial data yQ ~ ey^. 

Z°{x) = exp ( / yo{s)ds 



(15) 



Hence we see that it will not be possible to carry on explicit computations if we do choose y{t, 0) = H. 
Indeed, in that case, we would not know explicitly yx{t,0) (which is needed to compute the solution 
to system (HU). However, we are confident that this term is very small. Hence, we are going to go 
the other way around: we will choose our control explicitly in the Cole-Hopf domain and use it to 
compute our control «(■). Therefore, we are interested in the following heat system: 



SI 



Zt — eZxx — " 

Z(t,0) =1 

Zxit,l) =0 

Z{0,x) = Z°{x) 



Z 



on 
on 
on 
on 



(0,T) X (0,1) 
(0,T), 
(0,T), 
(0,1). 



(16) 
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If we go back to the Burgers' domain, this means that we somehow use the following boundary 
condition at x = 0: 

yAm = j^{fit,0)~H'). (17) 

We expect that the solution Z will converge towards H^{-), where H^{-) is the Cole-Hopf transform 
of the steady state /i^ : 

^ cosh (1^(1 -x)) ^ ^ 

= \,H ■ (18) 

cosh 

Indeed, we have the following lemma. 

Lemma 5. Let T > and Z^^ e H'^ (0,1) such that ^"(0) = 1 and Z°{1) = 0. Then system {I^ has 
a unique solution Z in the space L^({0,T); H^{0,1)) Ci H^({0,T); H^{0,1)). Moreover, there exists a 
constant C{Z'^) > depending only on \\Z'^\\fji such that: 

\\ZiT, .) - i/^(-)llHHo,i) < e-i/2C(Z°)e-^. (19) 

Proof. It is classical to show that system ([TO)) has a unique solution in the space ^^^((O, T);H^{0, 1)) n 
H^{{0, T); H^{0, 1)). One can even get more smoothness if needed. An efficient method is the semi- 
group method that one can find for instance in [35] . To compute the dynamics of system (|16p . we 
introduce the adequate Fourier basis of L^: 



}n{x) — a/2 sin iyjy'-^ + 2^ ""^'^ — 0- 



Hence /„(0) = /^(l) = 0. We will use the notation A„ — [n + Thus, /" = —^nfn- Let us give 
the following scalar products, which can easily be computed using integration by parts: 

{HVn) = 1^1^, (20) 

\(Z'IU)\ < ^(i + 5ll^1l„.)- (^i> 

In these equations (-I-) denotes the standard scalar product in L^{0, 1). Let us write Z ~ 1 + w. Hence 
w will satisfy w{t,0) = Wx{t,l) = 0. Easy computations lead to the following ordinary differential 
equations for the components of w on our Fourier basis: 

Wn{t) = ~e (^Xl + Wn{t) - ^(l|/n). 

It is easy to see that the fixed points for these ODEs are the expected coefficients (iJ^ — l|/n). We 
can solve these ODEs with our initial condition: 

Wn{t) = a„e-K^" + 5^)* + {H' - l|/„), 

where: 

an^{Z"\fn)-{H^U). 

Now we can estimate Z{T, •) — H^{-): 

ll^(T, •) - i/^(-)II^Ho,i) = E 

ri>0 

From the expression of a„, (|20p and (PTjl we get the easy bound: 

\2 „,2 ^ 1 « I II '70||2 



A>^<16+||Z"||^,(o^i), VneN. 
Thus, we get 

\\Z{T, .) - i/^(-)ll?^i(o,i) < (l6 + \\ZYhho,i)) e-"^ E 



n>0 
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Now we split the sum in two parts: n < N — [1/eJ and n> N. We get: 



|1Z(T,.)-F^(-)I1hmo,i) 




Hence, for e small enough, we have: 

||Z(T, .) - H^■)\\mio.l) ^ (J + l) + WZYmioA)) e""^- 

This concludes the proof of Lemma [5l □ 

Now we can prove Lemma |31 

Proof of Lemma^ Definition of the control: Using Lemma [SI we start by considering the 
solution Z e L2((0,T);i73(0, l))r\H^{{0,T); H^{0,1)) to system (UHl) with the initial data Since 
Zo{-) > 0, the usual strong maximum principle (see [38]) guarantees that Z{t,x) > 0. Thus, we can 
define: 

Hence y G X is a. solution to (O with initial data eya and boundary control v{t) — —2£Zx{t, 0). Since 
Z e L2(-(-o r).^3(o, 1)) ni?i((0,r);i/i(0,l)), we can show that its boundary trace Z^(i,0) belongs 
to H^/^{Q,T). Hence v e H^/^{0,T). 



Proof of an L°° bound on the solution: If e is small enough, then e||yo||oo < H. Moreover, 
we know that v £ H^^^{0,T). Hence, v e C"[0,T]. Assume that supjg t] ^ > Let Tg be a time 
such that v(Tq) — sup[Q j-^v > H. On the one hand, by the comparison principle from Lemma [TJ we 
know that: 

y<viTo) on (0,T)x(0,l). (23) 
On the other hand, we recall relation ([T7)) : 

y-(i>o)-^(f(i,o)-ff2)^ 

Hence, since v{To) > 0, we get yx{TQ,0) > 0. Thus, there exists a; > such that y{To,x) > v{To) = 
sup[o r] This is in contradiction with assertion (|23p . Hence, if e is small enough, «(•) < H and 

yiT, ■) < H. 



Derivation of the estimate at time T: Now we want to prove estimate (fT5)) from 
Lemma U) We want to use estimate ([TUl) from Lemma [S) We perform the following computation at 
time T and for any x G [0, 1]: 



2e 



Z 



= 2e 
< 2e 



Z{Z^-Hl)+Z^{H^ -Z) 



Hi 



ZH^ 



2e 



Z -H^ 



Thus, we get: 



||y(T, •) - /i^(-)|lL.(o 1) < (2e + MT, x sup — x ||Z(r, •) - H%-)\\^,^, 

^ ^ [0,1] ^ 

Now we use that \\y{T, ■)\\ac < H and sup[Q l/H^ < e+^Z'^". Hence , using also ((T9)) . 



|ly(r, .) - /i^(-)llL2(oa) ^ 7f + ^) C{Z'')e 



This estimate concludes the proof of Lemma 2] 



□ 



Remark 1. In Lemma^ we take an initial data yo G Ho{0, !)• This is a technical assumption that 
enables us to use stronger solutions. We will get rid of it later on, by letting the Burgers' equation 
smooth our real initial data which is only in L^{0, 1). 



10 



3.3 Second step: going back to the null state 



Once we have reached the steady state /i^, we wish to go back to the nuh state. This is done by applying 
a suitable negative interior control u. The control v will only be following the global movement. The 
intuitive idea is to apply some negative control u on [0,r] such that u{t)dt = -H. Thus, we hope 
to reach some state that is below and above a boundary residue — H. However, this last statement 
is only true up to some small function (small as T — > 0). The key will be to choose the duration 
T of this step small enough (with respect to e). 



H 




Focus 




X = 0.9 



x = l x = 0.97 



X = 1 



Figure 4: Numerical simulation of the push-down towards the null state and the creation of a boundary 
residue. The final state y(T, •) is almost above the residue fc^(-) — H. 



Lemma 6. Let e > and H > be given data. Assume that 2e < H . We consider the evolution of 
an initial data yi G L^{0, 1). For any T > 0, we consider the following controls for t G [0,r]; 



H 



u{t) = 

vit) = H + [ u{s)ds. 
Jo 

Then the associated solution y G X to system satisfies: 

y{T,-)-k'{-) + H>6{T,-), 
where S Cz X is the solution to some Burgers-like system given below and is such that: 



\\6{T,-)\\,.<e"'^^'^\\y,-k^ 



L2 



2H 



(e- 



T/2e 



- 1 



(24) 

(25) 

(26) 
(27) 



Proof. Let T > and consider the controls defined by ([Ml) and (P5|) . Let us consider the associated 
solution y E X to We compare y to the solution z E X to the following system: 



zt + zzx - ezxx = u{t) in (0, T) x (0, 1), 

z(i,0) =v(t) in (0,T), 

z(t,\) ^v{t)-H in (0,T), 

z{Q,x) ^yo{x) in (0,1). 



(28) 



The comparison principle from Lemma [T] tells us that y{T, ■) > z(T, ■). Now we want to derive precise 
estimates for the solution z E X. We write: 



z{t,x) — k^{x) + / u{s)ds + S{t,x) 



(29) 
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where S E X is thus the solution to the following system: 

' St-e6^, + k'6, + (^6 + J^u{s)ds^ {k' + 6)^ =0 in (0,r)x(0,l), 

6{t,0) =0 in (0,T), (30) 

d{t,l) =0 in (0,T), 

(5(0, x) =yi{x)-k%x) in (0,1). 

Note that it is convenient in this proof to use fc^ in order to get exact zero boundary conditions 
S{t,0) = 6{t,l) ~ 0. We multiply the evolution equation of (|5D|) by i5 and integrate by parts for 
X G [0, 1] to get a i^-energy estimate on S: 



ld_ 
2 dt 



1 fi 

2 , , / r2 



S'+e I 6% = -J ¥55^ - / [5+1 u{s)ds ) + 5) J 

= IJ^Hk'h- I [S+ I u{s)ds]S{k% 
1 



2 , 

Now we use definition ([TTt and the assumption 2e < H: 



S'ik%~ / iZ(s)ds / Sik')^. 



:,lloo - 2e - 2etanh(l)2 " e ' 



Moreover, u{s)ds < H. Hence, 



2dtJ^ - 2e io e \Jo J ^ ^ 

Let us denote E{t) — \\S{t, •)||l2. Hence, one has: 

E{t)<—E+—. (32) 

From Gronwall's lemma, we get: 

E{T) < {E{0) + 2H) - 2H. (33) 

This concludes the proof of Lemma [51 □ 

This is the end of the hyperbolic stage. We need to perform the reverse scaling of @ to go back 
to y (and not y). We have shown that we are above some boundary residue — H . Hence, we have 
to study the evolution of the following initial data: 

^%x) = i {h'{x) -H) = ^ (^tanh (^^{l - x)^ - 1^ . (34) 

One should be scared by the size of this boundary residue that we are left with. Lideed, its 
size grows like l/\/s. However it has the important feature that its typical wavelength is e. Hence, its 
spectral decomposition will mostly involve high frequencies that will decay rapidly during the passive 
stage thanks to smoothing effects of Burgers' equation. 



4 Passive stage and dissipation of the boundary layer 

The goal of this section is to prove the following estimate concerning the dissipation of the boundary 
residue created in the previous section. Indeed, although its L^-norm increases as e goes to zero, 
regularization effects of the Burgers equation will dissipate it in any positive time T. 

Lemma 7. Let T > be a fixed positive time. For any e > 0, let us consider (jj G X the solution to 
the following system: 

Pt + (f>^x -^xx =0 m (0, T) X (0, 1), 

(j){t,Q) ==0 in (0,T), 

(j){t,l) =0 in (0,T), 

(/)(0,a;) =<^^{x) in (0,1), 
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where ^'^(x) is the boundary residue defined by ([^p. Then for any d > 0, we have the estimate: 

MT,-)\\mo.i)^0,^o{e'-'). (35) 



X = 0.5 

I L 



X = 1 
I 



y{t,x) 



\ - " 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 



/ ' 
/ ' 

I 

/ 



<^^{x) 



Figure 5: Numerical simulation of the dissipation of the boundary residue At time t — 0^ the 

boundary residue was of size ||fE'^(-)||co = 100. 



4.1 Cole-Hopf transform 



Once again, we are going to use the Cole-Hopf transform to derive precise estimates. Therefore, let 
us introduce the following change of unknown for x G [0, 1] and t G [0,r]: 

1 



z{t, x) = cxp y — — J (j){t,s)ds 

This leads to the following heat system for the new unknown z: 

' zt-z,, =0 on (0,r) X (0,1), 

z,(t,0) =0 on (0,T), 

z,(t,l) =0 on (0,T), 

z{0,x) ^Z^{x) on (0,1), 

where the initial data is computed from the initial data <i>'^: 

1 



(36) 



exp 



2 Jo 

l + ef("-i) 

H 

1 + er- 



<^'{s)ds 



(37) 



An important remark is that < 0. Thus, by the comparison principle from Lemma[TJ < on 
[0, T] X [0, 1] and z > 1 on [0, T] x [0, 1]. The backwards Cole-Hopf transform will give us: 



m = -2- 



:{T) 



Hence, using the fact that z> 1, we will have the following estimate: 

|0(r,.)| <2|z,(T,.)|. 



(38) 
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All we have to do is to study the L'^-iiorm of Zx{T). To ease computations, let us introduce: 



w = (1 + e ' )zx, 



(39) 



such that w is the solution to: 



Wt - Wxx 

w{t,0) 
w{t,l) 
w(0, x) 










on (0,T)x(0,l), 

on (0,T), 

on (0,T), 

on (0, 1). 



4.2 Fourier series decomposition 

We use Fourier series to compute w{T, •). We will use the following Hilbert basis of made of the 
eigen- functions for the Laplace operator with Dirichlet boundary conditions on [0, 1]: 

Gnix) = V2sm{mrx) for n >1. 
Let us compute the decomposition of w(0, •) on this basis. We integrate by parts twice: 



(w(0, •)|e„ 



V2 — e ' [ sin(mTx)e " ^dx 
e Jo 



eV2 



sin(?i7ra;)e •= ^ 



- \/2r, 



cos(n7rx)e 







H 

\f2 enn 
H 1^^^ 



cos{mix)e " ^ 



-(^)'MO,.)|e„) 



Now we can estimate the size of u;(T, •) in L^(0, 1): 

\\AT,-)\\l. - 5](M0,.)|e„)-e-"^-^ 



ri>l 



< 



For a G R, the following easy inequality holds: 

. / 2 1 
— — < mm [ a ,- 

Hence we split the sum and cut at a level N{e): 

N-l 2 2 2 

iK^'-)iii. < 8E^+2E 

n=l k>0 



^-2{N+kY-K^T 



fc>0 

< +2 



We want to choose A^(e) — +00 such that e'^N^ — > 0. For instance, we can take N — [e^'^J , where 
77 > is small enough. For s small enough, we have: 

lk(r, < l^e^-^" + 4e-2^-^^^^^ = O (e^-^") . (40) 

Combining estimates (PH)) and (|38l) . and the definition we can easily deduce the estimate (|35p. 
This concludes the proof of Lemma [71 □ 
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4.3 Approximate controllability towards the null state 

First, let us prove the following technical lemma. Indeed, we have proven that the particular boundary 
layer $^ dissipates, but all we also want to know what would happen if we were very close to it. 

Lemma 8. Let us change the initial data from Lemma[7\ to ^^{x) + ^S^ . We assume: 

^'(x) + -S' < 0, (41) 
W{-)\\mo,i) = O,^o{s'). (42) 
Then, the conclusion of Lemma\^ still holds. 

Proof. We follow the same scheme than for the proof of Lemma [71 Hence, we start by taking the 
Cole-Hopf transform of the new initial data $^(x) + ^(5^. Therefore, after the Cole-Hopf transform 
we have the following initial data: 



From our previous computation ([37]) of Z^, we know that |(Z'^)j,| = 0{l/e). Hence, using condition 
we have: 

Z%x) ■ fexp r -l) = O.^oie). 

\ V ^£ Jo / / ffi(04) 

Let us use the fact that our heat system ([55)) is linear. Therefore, using the conclusion of Lemma [7| 
we have: 

|!z(T, ■)\\m(o,i) = O.^oie'-') + O,^o{e). 

Once again we apply the backwards Cole-Hopf transform. We use the fact that z > 1 (this comes 
from the comparison principle and the hypothesis (|4ip '). Hence, 

II0(t,-)IIl-^(o,i) <2|k(r,-)lki(o,i)- 

Thus, the conclusion ([55)) of Lemma [7| still holds with this new initial data. □ 

Now everything is ready for us to show the following small time approximate controllability result 
for system ([IJ. We have to combine the different estimates. 

Theorem 4. Let T,r > and yo E L^{0,1) be given data. Then there exists u,v G L°°{0,T) x 
H^I^{Q,T) such that the associated solution y Cz X to system on [0,T] satisfies: 



\\yiT,-)\\mo, 



< r. 



Proof. Take T,r > and yo G L^(0, 1) given data. Let us take a small e > and break down our 
time interval into four parts. We introduce Ti = T/3, T2 — Ti + e and T^ — T2 + e*. The first part 
[0, Ti] of length T/3 is designed to smooth the initial data. The second part [Ti, T2] of length e is the 
part where the settling of the boundary layer takes place. The third part [T2,T3] of length is the 
quick push down to zero. The fourth part [Ta.T] of length at least T/3 (when e is small enough) is 
the passive stage for the dissipation of the boundary layer. Let us give some details. 



Smoothing of the initial data: First, for t G [0,ri], we choose u{t) — v{t) — 0. The system 
evolves freely. Regularization effects of the Burgers' equation smooth our initial data yo e L^{0, 1). 
We have y{Ti, ■) G -ffd(0, 1). There are many ways to prove such a result. For instance, one can take 
the Cole-Hopf transform and use well-known regularization properties of the heat equation. 

Settling of the boundary layer: Next, for t G [ri,T2], we perform the scaling @. We want 
to apply Lemma [5| for a duration 1. Hence, let us choose some H such that H — 2 > 0. We take 
V G H^/^{0,T) the control from Lemma [H For t G [Ti,T2], we use: 

u{t) ^ 0, 

v[t) = -V 



e \ e 
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Controls used 


1 settling 


v{t) 

\ dissipation 




t = 


T 


L T2 ] 


1 T3 T 








1 u{t) 


smoothing 









time 



push-down 

Figure 6: Approximate null-controllability strategy. 
From Lemma 21 we know that: 

v{T2,-)--K{-) 



--3/2^-||(H-2 



)■ 



(43) 



Push-down towards zero: Then, still in the context of scaling (jH), we want to apply Lemma |S] 
during a very short duration . Hence, for t G [T2,T3], we choose the controls found in Lemma |5] 
(with a total time e^), and we scale them appropriately. That is to say: 

Combining and Lemma [SI we get that, at the end of this hyperbolic stage: 

> y{Ts, •) > + J<^(e^ , 
where (using estimate P^ ): 

\\5{e^-)\\^, + W-k^\\^=0,^o{e'). 



Dissipation of the boundary residue: Now we enter the passive stage. We choose v{t) = 
u{t) = for t e [TsjT]. Since e goes to zero, T ~ T3 > T/3. Hence we can apply Lemma|H]on a time 
interval independent of e. By using the comparison principle from Lemma [T] we can conclude that: 

||y(T,.)||i. =Os^o(e^-"), 

for any rj > 0. For instance, one can choose r] — ^. Then we choose e small enough to ensure that 
\\y{T, ■)\\j^2 < r. This concludes the proof of TheoremS) □ 

Remark 2. In the proof of Theorem^ we concatenate different controls found in different parts. This 
could be a problem for smoothness because we did not check compatibility conditions at the jointures. 
However, the proof provides a control v € H^^'^{0,T) and this doesn't require compatibility conditions. 
If one wants smooth controls, it is also possible. One can choose a smooth control close to our control 
for the approximate controllability, then end with a smooth control for the exact controllability. 

5 Parabolic stage and exact local controllability 

Theorem m takes care of the small time global approximate controllability towards the null state. To 
get Theorem[Tl we need to combine it with a small time local exact controllability result in the vicinity 
of the null state. We give in this section two different approaches for this type of result. 
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5.1 Fursikov and Imanuvilov's theorem 



The following theorem is due to Fursikov and Imanuvilov. Indeed, the techniques they expose in their 
book [U| can be applied to show the following result. However, the proof of this precise statement is 
not written, and one has to work to show that the control can be chosen to be smooth. 

Theorem 5. Let T > 0. There exists r > such that, for any initial data ya G i^(0, 1) satisfying: 



||yo||L2(o4) < 

there exists a control v G C^[0,r] such that the solution y E X to the system: 

yt 



yyx -yxx =0 

y{t,Q) =v{t) 

y(tA) -0 

y(0,a;) ^ yo{x) 



in (0,T)x(0,l), 

m (0,T), 

m (0,T), 

m (0,1), 



(44) 



(45) 



satisfies y{T, •) = 0. 



The full theorem is in fact more general since one obtains local exact controllability to the tra- 
jectories of system (05]). The proof relies on Carleman estimates for parabolic equations. It is an 
extension of a previous result with two boundary controls whose proof can be read in |24) . 



5.2 Using Cole-Hopf and a moments method 

In this section we give a proof of Theorem [5] (both for the sake of completeness and for avoiding 
Carleman estimates). It relies on the Cole-Hopf transform and a moments method introduced in |20] 
by Fattorini and Russell. 

Proof. Let T > 0. First, we consider the following heat system: 



z(i,0) 

^.(t,l) 

z(0,a;) 





a(t) 


zq{x) 



m 
in 
in 
in 



(o,r) X (0,1), 
(o,r), 

(0,T), 
(0,1). 



(46) 



This is typically a setting for which we can apply the moments method of Fattorini and Russel exposed 
in [20]. They prove this system is null controllable for any positive time by means of very smooth 
controls. Let us use some control a G C^[0,r]. They also prove that there exists some constant Ct 
such that the size of the control is bounded from above by Ct x H-^oHl^. Therefore, if zq is small 
enough in L^, one can steer it to zero with a control Q!(-) such that ||a(-)||oo < 1- 
Now we get back to our Burgers' system. For yo G L^(0, 1), let us choose: 



zq{x) ~ exp 



yo{s)ds 



Thus, if yo(') is small in L^{0, 1) then zo(-) too. If they are small enough, then we can steer zq to 
with a control such that ||a(-)||oo < 1- In that setting, we have z(-) > —1 thanks to the maximum 
principle for the heat equation. Hence, if we let y = —2zx/{l + z), we get a solution y G X to (pS)) 
such that y{T, •) = provided that condition (pi)) is satisfied for some r > depending only on T. □ 



6 Conclusion 

In our work, we want to underline two important ideas. The first one is the rigorous analysis of 
the hyperbolic limit system and of the adequate weak boundary conditions. These weak boundary 
conditions somehow describe the behavior of the boundary layer and what it will be able to do or not. 
The second idea is the dissipation of the boundary layer by the fluid system itself during the passive 
stage. Once a boundary layer is created, will the system be able to dissipate it in short time or not? 

These two ideas might be important for the analysis of more complex problems such as the Navier- 
Stokes Open problem [TJ For instance, one could try to see if the boundary layer appearing in [11] 
when trying to control the 2D Navier-Stokes system with Navier slip boundary conditions can be 
dissipated in small time by the system itself. 
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